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Abstract 

We show that the Q-box expansion of nuclear many-body physics can be applied in nuclear effective field 
theory with explicit pions and external sources. We establish the corresponding power counting and give 
an algorithm for the construction of a hermitean and energy-independent potential for arbitrary scattering 
processes on nucleons and nuclei to a given order in the chiral expansion. Various examples are discussed in 
some detail. 
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1 Introduction and summary 



Effective field theory methods applied to few-nucleon systems have become a powerful tool to establish a 
systematic and precise formulation of nuclear physics consistent with the symmetries of QCD, in particular its 
spontaneously and explicitly broken chiral symmetry. The field was initiated by a series of papers of Weinberg 
0121 El and has matured considerably. Recently, two-nucleon dynamics based on effective field theory (EFT) 
with explicit pions # 2 has been worked out at next-to-next-to- next-to-leading order (NNNLO), see |S] # 3 . 
Due to the appearance of shallow bound states a purely perturbative treatment like in chiral perturbation 
theory for systems with pions or with one nucleon is inappropriate. However, as stressed already by Weinberg, 
one can apply the power counting in the construction of the so-called effective potential. Through iteration 
of this effective potential one obtains the full physical scattering amplitude. In the case of nucleon-nucleon 
scattering, it was shown in that the power counting is violated if one considers diagrams containing two- 
nucleon intermediate states. From that, the simplest definition of the effective potential can be given as follows: 
the effective potential is the sum of all time-ordered diagrams which contain in each intermediate state at least 
one pion. No violation of the power counting appears on the level of the effective potential, thus it can be 
constructed perturbatively in terms of the usual small parameters g/A x and M v /A x , with q a generic external 
momentum, M n the pion mass and A x ~ I GeV the scale of chiral symmetry breaking. The extension of these 
arguments to pion production or the inclusion of the A(1232) is straightforward if one takes into account the 
new scales appearing in these cases (see e.g. [H]). We will not consider such extensions here specifically, although 
our formalism is general enough to include these. Once the effective potential is constructed to a given order in 
the chiral expansion, the full T-matrix generating the bound and scattering states can be calculated by solving 
the Lippmann-Schwinger equation, T — V e s + V c sGqT^ with Go the two-nucleon propagator. Even though 
this definition of the effective potential is very simple, one encounters some disadvantages. The potential is 
energy-dependent and not hermitean. This makes it difficult to apply it to scattering processes involving more 
than two nucleons. A solution to this problem was given by Epelbaum and collaborators who showed that 
the method of unitary transformations developed by Fukuda, Sawada and Taketani [H] and by Okubo [H] can 
be applied to generate an energy-independent and hermitean potential consistent with the power counting of 
chiral perturbation theory (the effective potential). While this method is quite powerful, it suffers from the 
disadvantage that for different processes one has to define different model spaces and each time to construct 
the operator A, which parameterizes the corresponding unitary transformation, and to deduce the resulting 
hermitean effective Hamiltonian. We propose here another scheme based on the Q-box expansion of Kuo and 
collaborators jlOl ITT] . The main results of this study are: 

i) We have shown that the Q-box expansion respects the chiral expansion: only a finite number of Q~boxes 
contribute to the effective potential at a given chiral order. 

ii) We have constructed an explicit algorithm that allows one to construct the effective potential for any 
process involving nucleons, pions and photons (or other external sources) to a given order in the chiral 
expansion, see section [SJ 

iii) We have shown for various examples that when one requires the on-shell condition to the asymptotic 
states, one recovers the expressions based on time-ordered perturbation theory, as it should be. 

iv) An explicit application of the method developed here is the calculation of the fourth-order corrections to 
the three-body contributions in neutral pion electroproduction off the deuteron, see |12j . 

The manuscript is organized as follows. In sectionElwe briefly review Weinberg's approach to the chiral dynamics 
of few-nucleon systems and the corresponding power counting. A general derivation of the hermitean effective 
potential is given in sec. [31 following closely ^3]. The Q-box expansion is discussed in sec. 0J These three 
sections do not contain new results but are needed to keep our presentation self-contained. Sec. [5] contains our 
central results. We show that the Q-box expansion is consistent with the chiral expansion and further give an 
explicit algorithm to construct the hermitean effective potential to a given chiral order. Some leading order 
considerations are discussed in sec. El in particular the role of the so-called Okubo corrections (wave function 
orthonormalization diagrams). Various technicalities are relegated to the appendices. 

# 2 We do not consider the so-called pionless EFT in what follows. 

* 3 Note that only in [S] a detailed analysis of the theoretical uncertainties at this order is given as it is required for any sensible 
EFT analysis. 
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2 A brief review of Weinberg's approach to chiral few-nucleon dy- 
namics 



This section contains a very brief introduction to chiral perturbation theory with two (or more) nucleons and 
serves mainly to fix the notation. As pointed out by Weinberg |3 El some years ago, for processes involving 
two or more nucleons the standard power counting of chiral perturbation theory is valid only for a subset of 
diagrams. This subset consists of the time-ordered diagrams which in every intermediate state contain at least 
one pion (or delta). The sum of all the diagrams from this subset is called the effective potential, denoted by 
V^ff ■ The whole transition matrix or nuclear wave function can be computed numerically in the standard way 
by solving the inhomogcneous or homogeneous Lippmann-Schwingcr equation with the given effective potential. 
For a more detailed description, see e.g. [T^| 

If we are interested in pion scattering processes on nucleons or nuclei, we can obviously proceed in a similar 
way and define an effective potential composed of diagrams, which contain in every intermediate state at least 
two pions (or deltas). In this way we get an effective potential matrix 

y _ f (TrNN\V cfi \TrNN) (i:NN\V cfi \NN) \ 
off ^ (NN\V cS \nNN) (NN\V cS \NN) )' U 

where NN denote here a state with two or more than two nucleons. The elastic pion scattering potential 
can be decomposed in two parts as has been shown in Fig. ^ the first part contains all diagrams from the 
effective potential with the pion acting as a spectator. Let us denote it by (TrNN\Vjg\TrNN). The second part 
contains all other possible diagrams. It is obvious that the first part is given by the usual effective potential 
Vjg N multiplied by the three dimensional delta function of the pion spectator. If we want to describe scattering 
processes on nuclei, it is useful to define the well-known Lovelace operator U^J, called here the irreducible kernel 
K: 

K = ^ + V ^ E-H - VeB + ie ^ (2) 
where E denotes the total energy of the system and the potentials V* s and VJ S are given by 

V$E = H Q + V eB -HZ ) V/ ff = H +V ef{ -Hl, (3) 

rt 



and Ho, Hq and Hq denote the free Hamilton operators of the particles in the intermediate state, of the initial 
and of the final asymptotic states, respectively. In the case of pion scattering, e.g., V* e and vjg are given by 
the effective potential without the pion spectator contribution: 

V* -Vt -V ( ^NN\V$\kNN) 0\ 

V eS - V eS - Veff - I Q q ) • (4) 

All diagrams from K are called irreducible and all others are called reducible. To calculate the transition 
amplitude for scattering processes on nuclei, one has to convolute the irreducible kernel K with the corresponding 
nuclear wave functions. Since in every intermediate state of the kernel K the chiral counting can be applied (no 
infrared singularities appear if the nucleon mass approaches infinity) , one can calculate it perturbatively. The 
simultaneous expansion both in small momenta (masses) of Goldstonc-bosons and in one over the nucleon mass 
can be performed. The 1/m Taylor series expansion within the potentials can be done, since no unitarity cuts 
can appear there (note that we always define the model space such that the potentials are real). The cuts, 
which are responsible for the imaginary part of the scattering amplitude, are all outside the potentials (the 
inside/outside refers to the structure of the kernel given in Eq. (0), where the trivial 1/m Taylor expansion 
does not make sense. However outside the potentials one can perform the 1/m expansion within the energy 
denominator, which is symbolically given in the following equation: 

1 = 1 f5) 

E-E' + ie {E - E')o + A/m + B/m 2 + He' w 

Here (E — E')q denotes the corresponding energy difference E — E 1 in the static limit. In some cases the nucleon 
recoil corrections, denoted here by A/m, are sizeable and can not be neglected, for more detailed discussion 
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Figure 1: A nuclear potential decomposed into the pion spectator and the interaction part. 



of this topic see e.g. |16| . The power counting for a given diagram from the irreducible kernel K or nuclear 
potential is given by [3] 

v = l-N-2C + 2L + ^V i b ii Ai = di + ^ni-2, (6) 

i 

where N, C, L, Vi, ck and n., are the number of nucleons, of connected pieces, of loops, of vertices of type i, of 
derivatives (or powers of the pion mass) at the vertex i and of nucleons at the vertex i, respectively. In order to 
calculate a scattering process up to a given chiral order v max we have to sum up all the diagrams with v < v max 
from the kernel K and convolute them with the chiral nuclear wave functions fixed by the corresponding nuclear 
potential calculated up to order v max + 3. The added 3 here guarantees that the operators V^ s and K are of 
the same chiral order # 4 . 

Although all that looks nice, one well known disadvantage appears: the effective potentials defined in this 
simple way are energy-dependent and non-hermitcan. This makes it difficult to apply the formalism to the 
processes, where more than two nucleons are involved. Epelbaoum et al. [7] suggested another construction of 
effective potential using the techniques of unitary transformation, firstly introduced by Okubo [S] many years 
ago. They showed explicitly, how to construct the hermitean energy-independent 7V-nucleon potential up to 
a given order in chiral perturbation theory. To be consistent with this picture we have to extend the given 
construction to the scattering processes on nucleons and nuclei. The main purpose of this manuscript is to give 
such an extension, such that both the nuclear wave functions on the one hand and the irreducible kernel on the 
other hand were constructed using hermitean energy-independent potentials in a unified way. 



3 Hermitean effective potential 

To derive the general effective potential, we proceed in the manner of Suzuki and Okamoto . For earlier and 
related work, see e.g. [T7| CS d E01 EH E3 E3 ■ Let the full Fock-space H of all physical states be described 
as a direct sum of a model subspace TLm and its complementary space Tin, 

H = H M ®H n - (7) 

Let P and Q be projection operators onto these two subspaces, 

P,Q:H -> H (8) 

such that Im[P] C Hm, Im[Q] C Ti-iii PQ = QP = and P + Q = 1, with 1 = Id the Identity operator 
Id : H — > H. Let A : H — > TC be an operator with the following property, A\Hm\ C TC-jz and A [H-r] = 0, such 
that it maps the states from the model space TIm to the states from the complementary space Tin- From this 
property it follows immediately 

A = QAP , (9) 
AQ = , (10) 
PA = , (11) 

and 

A 2 = . (12) 

# 4 Note that V e f is given by a nuclear potential multiplied by the pion spectator delta function, which decreases the chiral 
dimension by —3. 
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Now we define an operator X(n) by 

X(n) = (1 + A)(l + A^A + AA r ) n , (13) 

with n a real number. The inverse of X(n) is 

X-^n) = (1 + A^A + AA^)~ n (l - A). (14) 

The above relation can be verified by using the fact that (1 — A)(l + A) = 1. The stationary Schrodinger 
equation of the original problem H^> = E^> , where H = Hq + V is the Hamiltonian of the physical problem 
represented as a sum of the free and interacting parts Hq and V, can be transformed to 

X-\ri)HX{n)X- x (ri)yi = £A" _1 (n)*. (15) 

Suppose that the transformed operator X~ 1 (n)HX(n) fulfills for n — the decoupling equation 

QX- 1 (0)i7X(0)P = Q(H + [H, A] - AHA)P = 0. (16) 

This is exactly the equation of Okubo (Hj for the operator A. The above equation leads to the general decoupling 
equation 

QX^irfHX^P = (Q + AA i )- n QX- 1 (0)HX(p)P(P + A i A) n = 0, (17) 
where we have used the fact, that 

(1 + A^A + AA*) n P = (P + A^A) n , (18) 
Q(l + A^A + AA r ) n = (Q + AA^) n . (19) 

If we now start from the effective equation in the model space 

PX- l {n)HX{ri)P<£> = E$, $ e H M (20) 

we get (because of the decoupling equation) 

X _1 (n)i/X(n)$ = E®. (21) 

Multiplying this equation by X(n) from the left brings us to the original Schrodinger equation: 

= EV, (22) 

with = X(n)$. Let us introduce at this stage some important notation for the effective Hamiltonian and the 
potential: 

fleff(n) := PX-^HXWP = PH P + V eS (n) (23) 

To study the properties of V e g(n), we multiply the Eq. I|16|) from the left by A 1 and conjugate the last equation: 

PA^QHP + PA^QHQAP - PA^APHP - PA^APHQAP = 0, (24) 
PHQAP + PA^QHQAP - PHPA^AP - PA^QHPA^ AP = 0. (25) 

Eliminating the QHQ-tcrm from the above equations we get 

P(l + A ] )HP(P + A ] A) = (P + A^A)PH(1+A)P. (26) 

The effective Hamiltonian and its adjoint have the form 

H cS (n) = (P + A^A)- n (l-A)H(l+A)(P + A^A) n (27) 

= (P + A^A)- n PH(l + A)P(P + A^A) n , 

Hl s (n) = (P + + A*)H(1- A^)(P + A* A)~ n (28) 

= (P + A^A) n P{\ + A*)HP{P + A*A){P + A^A)-"- 1 

= (P + A^A) n+1 PH{l + A)(P + A^A)- n ^ 1 . 
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So we obtain the relation 

H cS (n) = fft s (- n -l). (29) 

From the above relation we see that H e f{(n) becomes hermitean when n = —1/2. 

For practical purposes two effective interactions are of special importance. The first one is i? o ff(0), which 
has a relatively simple form 

H eS (0) = PHP + PHQA = PH P + R, (30) 

with the non-hermitean potential, denoted by 

R = Vcfr(O) = PVP + PVQA. (31) 

Because of its simplicity this energy-independent non-hermitean potential has rich applications in nuclear 
physics. 

Another important effective Hamiltonian is H e g(— 1/2), which is hermitean. Its explicit form is 

H eS (-l/2) = (P + A^A) 1 / 2 H{P + A){P + A^A)- 1 / 2 = PH P + W, (32) 

with the hermitean potential, denoted by W = V c s(— 1/2). The hermitean theory is more complicated than the 
standard non-hermitean theory, but has an advantage that the wave functions, which come from the solution 
of the effective Schrodinger equation, are orthonormal. To see this, let us define an operator 

G = arctanh (A - A 1 ) . (33) 

It can be shown that the following relation is valid 20 : 

exp(G) = (1 + A - A+) (1 + A*A + AA^)~ 1/2 . (34) 
Using Eq. (jTHjl we immediately see that 

exp(G)P = (1 + A) (P + A^Ay 1 ' 2 = X(-1/2)P. (35) 

Let {iV'n)} denote a set of orthonormal eigenstates of the Hamiltonian H . The corresponding eigenstates {|</> n )} 
of effective Hamiltonian are related to the original one by 

|V„) = X(-l/2)|0n) = X{-l/2)P\4> n ) = exp(G)P|0„), (36) 

such that we get for the effective eigenstates 

\<j> n ) = exp(-G)|^ n ). (37) 
The orthonormality of the effective eigenstates follows now directly from the unitarity of exp(G): 

{<j>n\<i>m) = (ipn\ exp(G) exp(-G)|^> m ) = (4>„\ip m ) = S„. m - (38) 

Now let us discuss some interesting relations between hermitean and non-hermitean potentials. One can show 
the following general relation [T§| 

oo 

W = PVP+ J2 F{m,n)((A**A) m PVQA{A^A) n + h.c.) , (39) 

771,71—0 

with the function F(m, n) denned as coefficients of the Taylor series of 



r\ 

f(x,y) = -=L±?— = F(m,n)x m y n . (40) 



G 



Table 1: The values of F(m, n). 



(m, n) 


(0,0) 


(0,1) 


(1,0) 


(0,2) 


(1,1) 


(2,0) 


F(m, n) 


1/2 


-1/8 


1/8 


1/16 





-1/16 



From the property f(x, y) + f(y, x) = 1 we have the anti-symmetry relation 

F(m,ri) = -F(n,m), (41) 

unless m = n = 0. The values of F(m,n) for small m and n are given in Tabled Due to the anti-symmetry 
relation (I4H and from the definition of R in Eq. l|31|l , we can write W in another form as 

w = R+^W + ^ F{m,n){A^A) m {R-R)){A^A) n , (42) 

m+n>0 

where the summation runs over zero or all positive integers except m = n = 0. The expression of W in the above 
equation clarifies how the non-hermiticity term ~ (R — R>) contributes to the hermitean effective interaction 
W. The first two terms in the expansion of W are given by 

W = l-(R + R^ + l((A^A)(R-R^ + h.c.)-^-({AU) 2 (R-R r ) + h.c.)+... . (43) 
2 8 lb 

This formula will be useful for our later investigations. 

4 The Q-box expansion of the effective interaction 

We are now interested in the solution of the decoupling equation l|16fl for the operator A. Once we have 
constructed the operator A we are able to give an explicit form for the effective potentials R and W. For this 
reason we define a so called Q-box, which was first introduced by Kuo and collaborators JT1| in the study of the 
folded diagram expansion of an effective potential: 

Q(E) = PVP + PVQ — - 1 q H q QVP. (44) 

The physical meaning of this object is very familiar. It is nothing else but the sum of all time-ordered diagrams, 
which do not contain any model-space intermediate states. For example if one fixes the model-space as a space 
of all pure nucleon states, then every diagram from the Q-box has to contain at least one pion or delta in every 
intermediate state. Some lowest order diagrams of a Q-box in the case of nucleon-nucleon scattering are given 
in Figure We also want to introduce higher Q-boxes as 

C„( El ,..., E „ +1 ) . ^rPVQ (Ei _ QHQ) ,\ En+i _ QHQ) QVP. (45) 
Higher Q-boxes can be easily derived from the original Q-box through the relation (221 

n+l 

Q n {E x ,...,E n+l ) = J2 C k( E i>-> E n+i)Q( E k) (46) 
fc=i 

where 

fe ^ l ' ' n+1 {Eh — E x ) ■ ■ ■ (E k — E k _x){Eu — Ek + i) ■ ■ ■ (E k — E n+i ) ^ ^ 
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Figure 2: Some lowest order diagrams of a Q-box in the case of nucleon-nucleon scattering. 

If at least two of the energies mentioned above are equal, one has to understand Eq. I|46(l as a limit. To make 
this clear, let us look for example at Qi(E±, £^)-box, with the energies E\ = E^- 

+ p^^(_L___L_)q F p 

= lim M PVQ ^-L_ ( — _Q H Q^ }Q VF 



1 

(# 2 -~QHQf 



With these definitions one can express the potentials i? and W in terms of Q-boxes. Let us classify all 
contributions to W by a k number given by 

k (Qh ■ ■ ■ Qk n ) = k x + ■ ■ ■ + k n . (48) 

Denote by W n all Q-box contributions to W with k = n. The potential can be expressed as an infinite sum of 
those contributions: 

oo 

W = ^W n . (49) 
n=0 

Here we want to derive explicit expressions for Wo, W\ and Wi. For this reason let \a) be an eigenstate of 
PHqP with an energy E a from the model space, then the decoupling equation (|16[) can be written as 

Ala) - ^-W^-^W^- (50) 

For brevity let us now adopt a notation that whenever a symbol |/3)(/3| appears the summation over (3 has to 
be understood. Now we eliminate A from the potential R: 

R\a) = PVP\a) + PVQA\a) 

= Q(E a )\a) - PVQ Ea ^ QRQ A \P) (P\R\<*) 

= Q(E a )\a)+Q 1 (E a ,E M{m^+PVQ Ea _ 1 QHQE ^^ AMfrjiZjflffjiija) (51) 
= Q{E a )\a) + Q 1 (E a ,E )\(3)(0\R\a) +Q 2 (E a ,E l3 ,E J )\>y)(>y\R\0)(l3\R\a). 



8 



An iterative solution of the previous equation can be given by 

R\a) = Q(E a )\a)+Qi(E a ,E )\0){p\Q{E a )\a) 

+ Q 1 (E a ,Ep)\fl(p\Q 1 (E a ,EJ\7)('Y\Q(E a )\a) (52) 

+ Q2(E a , E p , E^\ 7 ){ 7 \Q(E )\^{/3\Q(E a )\a) +.... 

To derive an approximate equation for the operator A we implement the previous equation for R in Eq. l|5Uf) 
and get 

A|a) = ^^QVP\a) - ^-^A m p\Q(E M , (53) 

where we only take the first term of R , because we are not interested in the calculation of A, but of A' A (the 
neglected terms of R would lead to contributions with k > 2). The iterative solution of the above equation is 
given by 

A '°> = E-kmQ QVP ^ ~ E a - QHQ E f j \ HQ QV p mmE a )\*). (54) 
For the operator A* A we get 

(5\A<A\a) = -(5\Q 1 {E s ,E a )\a)-(S\Q 2 {E a ,E p ,E s )\p)(f3\Q(E a )\a) (55) 
- (S\Q(E s )\f3)(l3\Q 2 {E ai E ,E s )\a). 

If we now use the expressions for the operators R and A* A in Eq. I|43[l (note that the Q-box expansion of A A 
starts with a Qi-box such that all terms of Ea. (|42|l with m + n > 2 lead to contributions with k > 2.), we get 
the following explicit form of Wo, W\ and W 2 : 

(5\W Q \a) = ~{{5\Q(E a )\a) + {8\Q(E s )\a)\, (56) 
(S\Wi\a) = ^{5\Q x {E a ,E )\p)^\Q{E a )\a) + (5\Q(E s )\fl) (f3\Qi{Ep,E s )\a) 

+ -{-(SlQ^Es^EffM ((f3\Q(E a )\a) - (f3\Q(E p )\a) 

- ( K {5\Q{E s )\p)-{5\Q{E p )\p)Yp\Q 1 {Ep,E a )\ a )^, (57) 



2 v 

{Sl^iE^E^E^) ^1^)1/3) (0\Q(E a )\a) + (5|Q(f; 4 )|/3)(/3|Q(S /3 )|7)(7|Q 2 (f; 4 ,^,£; 7 )|a) 

^-{S\^i(E Si E p )\0) ((/3|Q 1 ( J B Q)J E; 7 )|7)(7|Q(^)|a) - (P\Q(E )\j)(i\Qi(E p , E^a)) 
(i8\G(Es)\0)(P\Gi(Es,Ef,)\>y) - {5\Qi{E 1 ,E p )\P)(f3\Q{E 1 )\ 1 )) ( 7 |Qi(£ 7 , E a )\a) 
{5\Q 2 {E 1 ,E p ,E s )\P)(f3\Q{E 1 )\ 1 ) ^\Q(E a )\a) - ( 7 |Q(£ 7 )|a)) 
(iS\Q(E s )\p) - (6\Q(E^m) (l3\Q 2 (E a ,E 1 ,E f3 )\ 1 )( 1 \Q(E a )\a) 
(SlQiEsmiPlQ^E^E^Eg)^) (( 7 |Q(E Q )|a) - ( 7 |Q(£ 7 )|a)) 

((6\Q(Es)\(3) - {5\Q(E P )\P)) {l3\Q{Ep)\^){^\Q 2 (E a ,E 11 Ep)\ot)^ 
^USlQxiEs^pmiPlQxiE^E^) ^\Q(E a )\a) - ( 7 |(§(B 7 )|a)) 

((6\Q(Es)\(3) - {5\Q(E P )\P)) (0\Q 1 (E^,E 1 )\^){'y\Q 1 (E 1 ,E a )\a)\. (58) 



Figure 3: Triangle nucleon-nucleon diagram. The ordinary time-ordered reducible diagram on the left-hand- 
side becomes an irreducible one by shifting the left (dot-dashed) cut. One gets an irreducible diagram on the 
right-hand-side with one insertion from a Qi-box. 

Until now we did not specify the model-space, so the above equation is quite general and can be applied to 
every problem in chiral perturbation theory with more than one nucleon. 

5 Hermitean potential up to a given order 

The main advantage of the Q-box expansion is that it respects the chiral expansion: only a finite number of 
Q-boxes contribute to the effective potential up to a given chiral order. 

To demonstrate this issue, we fix the model space as consisting of states with two or more nucleons. The 
states with at least one pion or delta are assumed to belong to the complementary space. Now we can give a 
diagrammatical interpretation of Eqs. i|5t)|) . (|57|l and l|58|l considering nucleon-nucleon scattering. The interpre- 
tation of Eq. (|56|l is already given in Fig. [3 To give an interpretation of the other two Eqs. (|57|) and l|58|) consider 
first the triangle diagram of Fig. [3J By shifting the NN cut to the NNir cut we get an irreducible diagram 
with one insertion from a Qi _ box such that this diagram contributes to the term (S\Q(Es)\(3) (f3\Qi (Ep, E$)\a). 
Even more all possible contributions to {8\Q(E$)\(3)(l3\Qi(Ep : E^)\a) are generated in the same way: every 
contribution of it can be visualized by a definite diagram with exactly one NN intermediate cut which has to 
be shifted to the complementary space intermediate states in every possible way. The same is valid for terms 
including higher Q-boxes: all contributions of terms with k = n can be visualized by definite diagrams with 
exactly n NN intermediate cuts, which have to be shifted to the complementary space intermediate states in 
every possible way. 

From this interpretation we see that the higher the K-number of a Q-box contribution to the potential W 
is, the higher is the number of NN cuts which have to be shifted and consequently the higher is the number 
of loops which increases the chiral order of a given contribution. For this reason there is only a finite number 
of Q-box contributions if the chiral order is fixed. Although we restricted here our considerations to nucleon- 
nucleon scattering, the last statement is quite general: we could choose other model and complementary spaces 
and show with the same arguments that the Q-box expansion of the effective potential W stops after a finite 
number of terms if the chiral order is fixed # 5 . 

Now we are at the point of showing how one can construct the effective potential W up to given chiral order 
diagrammatically. The algorithm is indeed very simple: 

• Draw all possible covariant diagrams (reducible and irreducible) , which contribute up to a given "naive" ^ 6 
chiral order. 

• Draw for every covariant diagram all possible time-orderings and classify them by numbers (/-c-numbers) 
of model space intermediate cuts. 

• In each time-ordered diagram with n — n shift all n model space intermediate cuts to the complementary 
space cuts in every possible way. Every diagram generated in this way represents a contribution to a 
definite Q-box structure with n — n. Summing up all these contributions to the energy propagator we get 
the irreducible energy denominator of a given time-ordered diagram. 

• The energy denominator of a given covariant diagram is then given by the sum of contributions from all 
possible time-orderings. 

* 5 A more formal proof for pion scattering processes off nucleons is presented in the appendix 

"naive" we mean here that we do not take into account the violation of the power counting by the reducible diagrams. 



to 



Figure 4: Triangle nucleon-nucleon diagram of chiral order v — 2. The two steps needed to determine the energy 
denominator are shown. For the discussion see the main text. 



• If K > and there are no possibilities to shift the model space intermediate cuts, then the corresponding 
diagram does not contribute to the effective potential W. 

The diagrams constructed according to this procedure obey the power counting, cf. Eq.©. To make the 
procedure clear let us discuss some examples: 

• The first one is the triangle diagram in nucleon-nucleon scattering shown in Fig. The model space in 
this case is generated by pure nucleonic states and the complementary space by states with at least one 
pion or delta. There are only two possible time-ordered diagrams with k = 1. The only possibility for 
shifting the nucleon-nucleon cut is to shift it to the 7riViV-state (see Fig. . All possible contributions to 
the potential W with k = 1 are given in Eq. Q57[). In this case the states |a), \j3) and \S) denote the two 
nucleon states \N3N4), \N$Nq) and \N1N2), respectively. Consequently 

E a =E 3 + E 4 , E = E 5 + E 6 , E S = E 1 +E 2 . (59) 

The energy denominators for the upper and lower time-ordered diagrams are given by 

Pu = -~ 1 - (60) 

2 E\ + E2 — E§ — £7 — £4 £/5 + Eq — — Ef — E4 



and 



2 Ei + E 2 — E 3 — £7 — E e E 5 + Eq — E 3 — E7 — E 6 ' ^ ^ 



respectively. The contribution weighted by 1/8 from Eq. (|57J) vanishes because the corresponding Q-box 
difference is of vertex type: 

(/3\Q(E a )\a) - (l3\Q(E p )\a) = (f3\V\a) - (f3\V\a) + ...., 

where (/3|V^|a) denotes the NN vertex. Ellipses denote here terms of higher order, which don't contribute 
to the considered diagram. The energy denominator P for the covariant triangle diagram is then given by 

P = P u + Pi. (62) 



• As the second example let us consider one-loop pion production diagram shown in Fig. The model 
space in this case is generated by pure nucleonic states and nucleonic states with one pion or delta. The 
states including two or more pions or deltas belong now to the complementary space. As we can see 
from Fig. [S] there are again two time-orderings. The difference to the previous case is that in the upper 
diagram there are no complementary space intermediate cuts. For this reason the model-space-cuts of the 
upper diagram can not be shifted and consequently this time-ordered diagram does not contribute to the 
potential matrix W . In the lower diagram there is a two-nucleon-model-space-cut which can be shifted 
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Figure 5: One-loop pion production diagram of chiral order v — 2. There are two possible time-orderings, which 
generate 1 diagram by shifting the NN cut to the cut with two pions. 



to the two-pion two-nucleon complementary space cut. The only contribution to the energy propagator 
comes from this diagram. To give the energy denominator we can use Eq. i|57|) , since the k number of the 
remaining diagram is equal one: 

P=-l I I . (63) 

2 Ei + £2 — £5 — £3 — Eg — Ej Eq + £7 — £5 — £3 — Eg — £7 



• As the third example let us consider pion production one-loop diagram shown in Fig. [B] In this case 
there are three time-orderings with k = 1 and one time-ordering with k = 2. By shifting the model 
space intermediate cuts we get 6 diagrams with k = 1 and one diagram with n = 2. To give the energy 
denominators of the first 6 diagrams with k = 1 we can use Eq. (15 711 : 



2 E\ + E-2 — E\ — Eq — Ej — £4 Ei + E2 — E5 — Eg — Ej — E4 
1 



Eg + Ef + E4 — _£/5 — Eg — Ej — E4 

1 1 1 

8 Eg + Ei + £4 — Er 3 — Eg — Ei — E4 E§ + E% + E4 — E§ — Eg — Ej — £4 
1 1 



Ei + Ei — Ei — Eq ~ Ei — E4 Eg, + Ei + £4 — £4 — Eq — £7 — £4 
P 2 = -i I \ 

2 £5 + £3 + £4 — Ei — Eq — £7 — £4 Eg + £7 + £4 — Ei — Eq — Ej — £4 

1 



£5 + £3 + Ei — Eb — Eg — £7 — £4 

1 1 1 

8 E\ -\- E2 — Ei — Eq — Ej — E4 Eg + Ej + _£/4 — Ei — Eq — Ej — E^ 

1 1 



£5 + £3 + £4 — £5 — Eg — £7 — £4 £g + £7 + £4 — £5 — Eg — £7 — £4 

1 1 1 

2 £1 + £2 — £5 — Eg — Eq — £2 Ei + £2 — £5 — Eg — Ej — Ei 

1 



Eg + Eq + E 2 — £5 — £9 — £7 — £4 



2 Ei + £2 — £5 — Eg — Eq — £2 £8 + Eq + £2 — £5 — £9 — Eq — £2 
1 



•^8 + Eq + £2 — £5 — £9 — £7 — £4 



(64) 
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Figure 6: One- loop pion production diagram of chiral order v = 2 (in standard Weinberg counting). There are 
four different time-orderings, which generate 7 different diagrams by shifting the ttNN intermediate cuts to the 
cuts with two or more pions. 



p» = 4- 



1 



2 Ei + E2 — E$ — Eg — Eq — Ei Eg + Eq + Ei — E§ — Eg — Eq — E2 
1 



Pa = 



Eg + Eq + E2 — E5 — E3 — Ei — E% — E% 
1 1 



1 



2 E\ + E2 — E5 — Eg — Eq — E2 Ei + E2 — £5 ~ E3 — Ej — Eq — E2 
1 



Es + Eq + E2 — E$ — E3 — Ei — E% — E2 



Pi 



1 



1 



1 



2 Ei + E2 — E5 — Eg — Ej — £"4 Eg + Eq + E2 — £5 — Eg — Ej — E^ 
1 



£•8 + £7 + £4 — £5 — Eg — £7 — £4 

The whole energy denominator is then given by a sum of all seven contributions: 

P = Pi + ■ ■ ■ + P7. 



(65) 



(66) 



To give the expression for P7 we used Eq. (|58|l for the n = 2 contribution to W . Note that there are no 
contributions to P7 weighted by 1/8 and 1/16, since the corresponding Q-box differences are of vertex 
type. 

From the third example we see that the expressions for the energy denominators can be rather lengthy. For this 
reason it might be nice to have an independent check of them. This indeed can be easily done: If one requires 
the on-shell condition (the initial and final energy of a given scattering process are required to be equal), then 
every time-ordered diagram (reducible or irreducible) can be written in terms of hermitean potentials, including 
some iterations from the Lippman-Schwinger equation. This is illustrated in Fig. [7] for the one-loop reducible 
triangle nucleon-nucleon scattering diagram. The corresponding potentials are denoted by a gray shaded area. 
The energy denominators of both diagrams, the first one (which comes from the first iteration in the Lippman- 
Schwinger equation) composed of two contributions and the second one composed of one contribution from the 
hermitean potential, are given by 



Ti = 



1 



1 



1 



£ — £5 — Eg 2 V £5 + £6 — £5 — £7 — £4 



£5 - £7 - £ 4 



(67) 
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Figure 7: One-loop triangle nucleon-nucleon time-ordered diagram written in terms of hermitean potentials 
denoted by gray shaded area. 



and 

2 E — E§ — Ej — E4 E§ + Eq — E5 — Ei — E4 ' ^ ^ 

respectively. Here we use the on shell condition E = E\ +E2 = E3 + E4. If we now sum up the two contributions 
Ti and T2 we get the familiar energy denominator of the corresponding time-ordered diagram^ 7 : 

Tl+T 2 = p W F~F w 1 W p - ' ( 69 ) 

In the same way one should be able to reconstruct every possible time-ordered diagram. 



6 Some leading order considerations 

Let us briefly discuss the differences between some hermitean and non-hermitean potentials in the case of pion 
scattering on few-nucleon systems. We denote the corresponding potentials here by W and W, respectively. 
Let us now consider some diagrams, which contribute to the potentials W or W 1 at leading order in the 
1/m expansion, namely in the static limit m — > 00. If we assume that the incoming and outgoing pions of 
the potentials W or W' are not involved in further interactions in the kernel, then (because of the on-shell 
condition) their energies have to be equal in the static limit. In this case contributions to the potential W 1 can 
be identified with the k = contributions to W. All other contributions to W with k > are called Okubo-type 
corrections (these corrections are also called wave function orthonormalization diagrams |24|). In some cases 



,-1 



,-' q 

— i 



+ 



+ 



Figure 8: Time-ordered diagram with two Okubo-type corrections for pion scattering off two nucleons. 

the Okubo-type corrections cancel some time-ordered diagrams from W' . Consider e.g. the pion scattering 
diagrams shown in Fig. [SJ The first time-ordered diagram contributes obviously to W and W. The two other 
diagrams are Okubo-type corrections. The energy denominators of these three diagrams are given by 

Pl = ^' P2 = ^^' P3 = ~^' (?0) 

such that their sum cancels. Here E q ii denotes the energy of the rescattered pion. The same happens if we 
consider the pion-scattering diagrams on three nucleons, shown in Fig. |5J The Okubo-corrections cancel again 
the time-ordered diagram. Like in the case of three nucleon scattering ^ 1251 1241 [7] it is straightforward to see 
that the contributions to W in the case of pion scattering on three nucleons, shown in Fig. 111)1 vanish. 



* 7 Note that this calculation serves only to check the correctness of energy denominators in the effective potential. A strictly 
perturbativc approach makes no sense in the case of nucleon-nucleon scattering. 
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Figure 9: Pion scattering on three nucleons: time-ordered diagram with two Okubo-type corrections. The 
ellipses denote the other time-ordered diagram with the corresponding Okubo-type correction. 




Figure 10: Two diagrams that contribute to the hermitean potential in the case of pion scattering off three 
nucleons. 
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A The degenerate case 

The aim of the three following appendices is to give the formal proof that the Q-box contributions with k > 2 
(see Eq. (14811 ) lead in the case of pure nucleonic scattering to operators of chiral orders v > 10 — 3iV and in the 
case of pion scattering off nucleons to operators of chiral orders v > 7 — 37V. Such kind of operators contribute 
in the case of pure nucleonic scattering to the orders higher than NNNLO. In the case of pion scattering on 
nucleons they lead to the kernel contributions of order higher than 0{q ) in the conventional counting, see e.g. 
2(i for a detailed discussion. 

Let us assume for the moment that all states from the model-space have the same energy E. In this case we 
can express the potentials R and W in terms of Q-boxes in a compact way. The operator PHqP can now be 
replaced by EP and the higher Q-boxes by 



1 d n Q 
~\~dE<- 



n 



(A.l) 



I \ "+ 1 



= { - l)nPVQ \j^QHQ ] QVR 



The decoupling equation (|16fl for the operator A becomes 



A = - QVP - AR, (A.2) 

E - QHQ E-QHQ ' v ' 



with R as given in the Eq. I|31fl 

R = PVP + PVQA. 
Eliminating A from both equations above and using the definition of Q n we obtain 



R = J2Q n R n , (A.3) 



n=0 
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where we have used the notation Qo = Q for ft = 0. This equation can be solved exactly as fT5) 

oc 

^ = E E f(m 1 ,...,m k )Q mi ---Q„ lk . (A.4) 



k=l mi,...,m k 

The /(mi, ...,rrik) are coefficients defined by 



<W,o for = 1 

/(mi,...,m fc ) = <( <W,o<W,i for A: = 2 (a.5) 



where is the step function 

nl \ f 1 for X > , k n\ 

9{X) = { for x < 0. (A ' 6) 

Let us introduce a convenient notation for representing the series of R. We define a symbol [m]k, which denotes 
a set of k elements which are zero or a positive integer 

[m]k = (mi,...,m fc ). (A. 7) 

We further define 

Q[m]k = Q mi ---Qm k - (A. 8) 
Using these simple notations, we may write R as 

oc 

R = EE/(Hic)QH^. (A.9) 

fc=l [m] k 

The coefficient f([m]k) has an interesting property: 

/(Mki)---/(MfcJ = /(™,Hfci,-.H*: n ). ( A - 10 ) 

such that the operator iZ n can be written as 

oo 

Rn = EE/k^wOh^ (A.11) 

fc=0 [m] fc 

where we use the notation 

f(n,[m] ) = f(n). (A.12) 
Now we can solve Eq. I|A.2(I for A exactly: iterating Eq. (|A.2J) we have 

A = £(-!)" (^3^) QVPR n (A.13) 

oo oo / 1 \ n "'" 1 

= E E E (-!)"/(", H fc ) L-^— QypQ[m] fe . (A.i4) 

n=0fc=0[m] fc ^ H) 

Multiplication with gives 

oo oo oc oo 

A ^ A = -EEEEEE/^H^/^wo^H^p+.+iQWi. (A.15) 

p=0 g=0 fc=0 Z=0 [m]h Hi 
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Therefore, we can write A' A as a sum of definite combinations of Q-boxes. We see from the above equation that 
the coefficients in front of the Q-box combinations in the sum have the value or 1. They serve to exclude all 
impossible Q-box combinations. From the definition of the function / in Eq. (|A.5J| we see, that the coefficients 
in the sum of Eq. (TO^ are zero unless 

p + 77ii + ■ ■ ■ + nih — k and q + ni + • • • + n\ = I (A. 16) 

If we write A^A as 

oo 

AU = -^^%]r)0[m] r , (A.17) 

r=l [m] r 

we obtain that the coefficients /i([m] r ) are zero unless mi + ■ • • + m r — r. This and the similar property of the 
coefficients f([m\ r ), namely they are zero unless toi + • • • + m r _i = r — 1 and m r = 0, will play an important 
role for the later power counting arguments. 



B Power counting for Q-boxes 

Let us derive the power counting for the Q(E)-box following Weinberg's arguments. A general time-ordered 
diagram of the Q-box has the following form: 

where L is the number of loops, R is the number of energy-denominators, di is the dimension of a vertex i and 
Vi counts how often the vertex of type i appears in the diagram. Under / p , E p , I n and E n we understand the 
number of inner pion, external pion, inner nucleon and external nucleon- lines. For the chiral order v one gets 

v = 3L-R-I p + ^2Vidi-3(C + D p + D n -l), (B.18) 

i 

where C, D p and D n denote the number of connected, disconnected pion (pion spectators) and disconnected 
nucleon (nucleon spectators) pieces, respectively. Each connected and disconnected piece brings a delta function, 
which carries dimension —3. The overall delta function is not taken into account, for this reason (C+D p +D n — 1) 
appears in the above equation. The phase-factors l/^/2E 7V of external pions are also not taken into account. 
With the help of the topological identities 

J2ViUi = 2I„+E n -2D n , (B.19) 

i 

J^ViPi = 2Jp + E p — 2D p , (B.20) 

i 

and 

L = C + L p + I n -J2 V *> ( R21 ) 

i 

R = J^Vi-l, (B.22) 

i 

we get the master formula^ 8 

v = A-3N - Dp- Ep + ^Vim, (B.23) 

i 

# 8 For Dp = one can easily verify that the equations IB. 231 and JS} are identical. 
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where 

3 

n t = di + -th + Pi - 4, (B.24) 

and rii denotes the number of micleons in vertex i, which can only be even, and pi counts the number of pion 
fields associated to the vertex i. 

In the case of the nucleon-nucleon interaction there are no external pions so that D p = E p = and we get a 
simplified formula for the chiral order 

v = 4-3iV + ^V*Ki. (B.25) 

i 

Because of chiral symmetry, has to be a positive number, thus the order of every diagram is bounded from 
below. The lowest possible order for a Q-box is v min = 6 — 3N and is represented through the lowest contact 
term from C^) N or the pion exchange diagram with the vertices from We understand from now on the 

Q-box as organized in a chiral expansion, 

oo 

Q(E) = J2Q l6 ~ 3N+i) ( E )> ( R26 ) 
»=o 

where the upper index denotes the chiral order. Analogously we can write the higher Q-boxes as 

oo 

Q„(£i,...,£„+i) = Y,Qn- 3N - n+l) (Eu-,E n+1 ), (B.27) 



i=0 



with 



n+1 



Q£- 3N - n+l) (Ei,-,E n+1 ) = J2 C k( E ^-^ E n+i)Q (6 - 3N+i) (E k ), (B.28) 

k=l 

where Ck(Ei, E n +i) is defined in Eq. I|47|l . The following combination of Q-boxes 

6 = Q^~ 3N ~ n+i) (E 1A , E ltn+1 )Q^- 3N - m+j \E 2 ,i, ...,E 2 , m+ i) 
= Qg- 3N - n + i \E 1A ,...,E w )\a)(a\Q%-^^ 

leads to the chiral order v{0): 

v{6) = 6-3N -n + i + 6-3N -m + j + 3(N - 1) 

= 9 - 3N + i + j - n - m, (B.29) 

where the factor 3(N — 1) comes from the integration over all relative coordinates of the intermediate state 
|a)(o;|. For a general combination of Q-boxes 

O = Qnx 3N (E%,1, Ei >rll+ x) ■ ■■ Q^ m 3N " m+Im ^(-E m ,l, — , ^m,n m + l) 

we get 

v{6) = to(6 - - (m + ... + Tim) + (h + - + i m ) + (m - l)3(iV - 1) 

= 3to + 3 - 3N - (m + ... +n m ) + («i + ... + i m ). (B.30) 

Now we are ready to proof that one can neglect higher Q-boxes for the case of pure nucleonic scattering. 
At the first sight it is not clear, why the Q-box contributions with increasing k number should increase the 
chiral order. Conversely, we see from the last equation that the order v(0) decreases with increasing k = 
rii + • • • + n m . To clarify this we first point out that the steps for the derivation of the operators R and A> A 
in the degenerate and in the non-degenerate case were the same. The coefficients in front of the Q-boxes are 
the same, the difference is only in the energy-structure of the Q-boxes. From the degenerate case we learned 
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that the coefficients f(n\, n m ) and h(ni, n m ) of the operators R and A' A restrict the numbers m, n m , 
namely f(n\, n m ) ^ only if the condition m + ... + n m — rn — 1 is valid and h(n\, n m ) ^ only if the 
condition n\ + ... + = to is valid. As a result we get for the operator R that 

ri\ + ... + n m = m — 1 (B.31) 

and from Eq. (|B.30(1 

j/(0) = 3m + 3-3N -m + l + (i 1 + ... + i m ) 
= 2m + 4-3N+(h + ... + i m ). 

The NNNLO will be reached when v = 10 — 3N. One obtains from the above equation that m < 3 and 
k = rii + ... + n m < 2. All the other possibilities are beyond NNNLO. For the operator A' A we get 

m + ... + n m = m. 

However, we have to remember that the operator A* A appears in the hermitean potential W only in combination 
with the operator R, which can be directly seen from Eq. Q42JI. For this reason a Q-box combination Q ni ■ ■ ■ Q nm 
from A' A has to be multiplied by the Q-box combination Qi ± ■ ■ ■ Qi k from R. So for the whole operator we get 
ni + ••• + n m = m and l\ + ... + Ik = k — 1 and altogether 

k = ni + ■•• + n m + h + ■•• + h = Tin + k — 1. (B.32) 

From the above considerations we see, that the condition m + k < 3 has to be valid. With similar arguments 
one can show that the operators of the form 

(AU) m R(A^A) n (B.33) 
with to + n > 2 contribute to orders higher than 10 — 3A^. 



C Power counting for pion scattering on nucleons 

In this section we want to consider the scattering of pions off nucleons, so we define our model-space as a space, 
which contains the states with only nucleons and nucleon states with at least one pion or delta. The states with 
more than two pions are in the complementary space. The power counting of Weinberg given in Eq. (|B.23(1 
remains the same but D p and E p do not vanish any more, because external pions can appear. 

We have again to prove the assumption that one can neglect the higher Q-boxes for this kind of processes. 
For this reason let us introduce a convenient notation: we write the projection operator P as P = P° + P , 
where P° and P 1 project to the pure nucleonic states and nucleonic states with one pion. Utilizing this notation 
one can write a Q-box as a two-by-two matrix: 



P 1 Q(E)P 1 P 1 Q(E)P° 
P°Q{E)P 1 P°Q(E)P° 



(C.34) 



The multiplication of Q-boxes is nothing but a matrix multiplication. At this stage there appears one problem: 
in the power counting Eq. © we did not count the phase-space factors of external pions. However, when we 
construct a Q-box we a priori do not know whether the outgoing pions of the box are internal or external. 
Indeed the assumption that all outgoing pion lines of a Q _ box are external results in a contradiction, because 
by multiplication of two Q-boxes external pion lines can become internal. For this reason let us introduce some 
kind of internal power counting, where we count all phase-space factors of pions, and denote the internal chiral 
order as v. The number of phase-space factors of external pions is given by E p — 2D p . We get the internal order 
by simply subtracting the factor [E p /2 — D p ) from the original order: 

v = v- -^-+D p = 4-3N- ^E p + ViKi (C.35) 



19 



Figure 11: P 1 Q(E)P 1 -boxes to leading order. 



Figure 12: P 1 Q(E)P° and P°Q( J B)P 1 -boxes to leading order. 

Every matrix- element has its definite minimal original and internal order. The minimal internal order of 
P 1 Q(E)P 1 is i> min = 3 — 3-/V and is represented by the diagrams of Figure ITTI The lowest internal orders of 
elements P°Q(E)P 1 and P 1 Q(E)P° is £ min = 3 — 3N +1/2 and is represented by the diagrams of Figure IT21 
Finally the lowest internal order of elements P°Q(E)P° is given by v XBiB = 6 — 3iV and is represented through 

the contact interaction. All the vertices in the above mentioned diagrams stem from C^n and £nn- To have 
an overview let us define an order matrix as 



. / P l Q{E)P 1 P 1 Q(E)P° \ j D^Q^P 1 ) v(P 1 Q(E)P°) 

V I P°Q(E)P 1 P°Q(E)P Q ) ~ y 9(P°Q(E)P 1 ) 9(P°Q(E)P°) 

With this notation we have for the minimal internal chiral order of Q-boxes 

. / P 1 Q(E)P 1 P 1 Q(E)P° \ _ ( 3-3iV 3-37V+± 
Vl ™ \ P°Q(E)P 1 P°Q(E)P° I ~ \3-3iV+i 6 - 3N 



(CM) 



(C.37) 



and for higher Q-boxes 



O = 



P l Q n {E u ...,E n+1 )P l P 1 Q n (E u ...,E n+1 )P° 
P°Q n (E 1 ,...,E n+1 )P 1 P°Q n (E 1 ,...,E n+1 )P a 



- / 3 - 37V - n 5-3N+±-n\ 

V - {0) = { 5 -3N + l-n 8-3N-n J' (C - 38) 

where n > 1. The maximal order, which corresponds to order q 4 in pion scattering off nucleons (standard 
single-nucleon counting) is given by v msx — 7 — 3N. The following matrix of Q-boxes 

/ plQ(3-™-n + i) {Eu ^ En+i) pl plQ(*-™^-^\ El> . 
I pOtf*-™H^+i) (Eli ^ E )P 1 P 0Q(?-™-^\ El ,... 
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can also be represented as 



P iiQQ-^*'W'^- 3N - n +i\ El ,...,E n+1 )P i >. 
For the order of m multiplied Q-boxes 

o = p <i ^ B r i(il+<a)+4ii - <2 " 3JV " ni+J ' l) (^i,i,-».^i,n 1+ i)i ,<a 

u i m A( 5 ~i ( - i ^+ i ^+^+ S ^.i m +i- 3N - n ^+3^')fj P T? \D»m + l 

we get 

v{6) = m(5-3JV)-(ni + ... + n rn ) + 0'i + -+jm) + (m-l)3iV (C.39) 

^ m m 

- +Wi) + Xi^.iMi -3(i2 + - + im)-3^5 ifc , . 

fe=l fe=2 

The last sum of Kroneckcr-symbols comes from the integration over the intermediate states. The integration 
over the relative coordinates brings in the case of a pure nucleonic intermediate state 3 (AT — 1) and in the case 
of a state with one pion 3N extra dimensions in the power counting. For the construction of the operator R we 
can use the following relations 

m-1 = m + .. + .n m , 

m 

m-1 = j 2 + ... + irn + ^ Sj k fl, 

fc=2 

where the fact that i 2 , ■■■,i m G {0> 1} has been used, and we get for the internal order 

( m 3 3 \ 

i>(6) = m-3N + 4+[^5 ik!ik+1 --i 1 --i m+1 \+(j 1 + ...+j m ) (C.40) 



2 2 

v.fc=l 



> m-3N + l, 



such that the relation m < 6 has to hold. However, the value of m can be further decreased. The above 
inequality becomes an equality, if i\ = i m +i = 1 and the sum of the Kronecker-symbols is zero. Let us discuss 
first this elastic case. For every k = 1, ...,m the condition if. ^ ik+i is valid. In this case at least one pure 
nucleonic intermediate state has to appear. It follows that D p — and the original order of this combination of 
Q-boxes is v = v + E p /2 = m — 3N + 2 leading to the condition m < 5. In the inelastic case i\ — or z m +i = 
the original order is given by v > m — 3N + 3, such that even the condition m < 4 has to be valid. Further 
decreasing of m will be obtained in the following manner: let us assume the equality m — 5. In the elastic case 
H = im+i = 1 the condition Y^k=i $ik,ik+i — ^ nas *° be valid. For m = 5 only the equality is possible, so that 
there appears at least one pure nucleonic intermediate state and the original order increases by one, such that 
v = 8 — 3./V and is beyond the fourth order. In the inelastic case the condition m = 5 is excluded and we get 
m < 4. Let us assume that m = 4. In the elastic case we get the condition Y^k=i&ik,ik+i — In any case 
one pure intermediate nucleonic state has to appear, otherwise we would get X^feLi ^i k Sk+i — 4- This restricts 
the above estimation to J2T=i ^*k>**+i — 1- F° r 171 — 4 on ly the equality YHc=i $ik,ik+i = is possible with the 
Q-box combination given by P 1 Q ni P Qn 2 P 1 Qn 3 P Qn 4 P 1 ■ However, for m — 4 at least one higher Q-box has 
to appear. For the above combination of Q-boxes one gets an increase of the original order by two, such that 
v = 8 — 3iV. In the inelastic case we get the condition J2T=i ^**>**+i = ^ wn i cn is impossible for m = 4. This 
proves that the Q-box contributions to R with n > 2 arc of higher order than 7 — 3A^, that is they start to 
contribute at 0(q 5 ) in the standard single-nucleon counting. The proof that the operators 

(A^A) m R(A^A) n (C.41) 

are of higher order than 7 — 37V if m + n > 2 can be given in the similar manner (the same holds for the k > 2 
contributions if m + n < 2). 
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